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Abstract —In this paper, an energy harvesting dual-hop relay¬ 
ing system without/with the presence of co-channel interference 
(CCI) is investigated. Specifically, the energy constrained mutli- 
antenna relay node is powered by either the information signal 
of the source or via the signal receiving from both the source and 
interferes In particular, we first study the outage probability and 
ergodic capacity of an interference free system, and then extend 
the analysis to an interfering environment. To exploit the benefit 
of multiple antennas, three different linear processing schemes 
are investigated, namely, I) Maximum ratio combining/maximal 
ratio transmission (MRC/MRT), 2) Zero-forcing/MRT (ZF/MRT) 
and 3) Minimum mean-square error/MRT (MMSE/MRT). For 
all schemes, both the systems outage probability and ergodic 
capacity are studied, and the achievable diversity order is 
also presented. In addition, the optimal power splitting ratio 
minimizing the outage probability is characterized. Our results 
show that the implementation of multiple antennas increases 
the energy harvesting capability, hence, significantly improves 
the systems performance. Moreover, it is demonstrated that the 
CCI could be potentially exploited to substantially boost the 
performance, while the choice of a linear processing scheme plays 
a critical role in determining how much gain could be extracted 
from the CCI. 

Index Terms —Dual-hop relay channel, wireless power trans¬ 
fer, co-channel Interference, linear multiple-antenna processing, 
performance analysis. 


I. Introduction 

Energy harvesting technique, as an emerging solution for 
prolonging the lifetime of the energy constrained wireless 
devices, has gained significant interests in recent years. The 
conventional energy harvesting techniques rely on the external 
natural resources, such as solar power, wind energy or ther¬ 
moelectric effects [1-3]. However, due to the randomness and 
intermittent property of external natural resources, communi¬ 
cation systems employing the conventional energy harvesting 
technique can not guarantee the delivery of reliable and 
uninterrupted communication services. Recently, the wireless 
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energy transfer technique, first demonstrated by Nikola Tesla, 
has rekindled its interest in the context of energy harvesting 
communication systems where radio-frequency (RF) signals 
are used as an energy source [4-8]. Since RF signals can be 
under control, it is much more reliable than external natural 
resources, hence, wireless energy harvesting using RF signals 
is a promising technique to power communication devices [9]. 

Since RF signals are capable of carrying both the infor¬ 
mation and energy, a new research area, namely simultane¬ 
ous wireless information and power transfer (SWIFT), has 
recently emerged. The seminal works [6,7] have characterized 
the fundamental tradeoff between the harvested energy and 
information capacity. Nevertheless, it was assumed in [6] that 
the receiver can decode the information and harvest energy 
from the same signal simultaneously, which is unfortunately 
impossible due to practical circuit limitations. To this end, 
the work in [10] proposed two practical receiver architec¬ 
tures, namely, “time-switching”, where the receiver switches 
between decoding information and harvesting energy, and 
“power-splitting", where the receiver splits the signal into two 
streams, one for information decoding and the other for energy 
harvesting. Since then, a number of works have appeared in the 
literature investigating different aspects of simultaneous infor¬ 
mation and energy transfer with practical receivers [11,12]. 
Specifically, in [11], an opportunistic RF energy harvesting 
scheme for single-input-single-output systems with co-channel 
interference (CCI) was investigated, where it was shown that 
the CCI can be identified as a potential energy source. Eater 
on, an improved receiver, i.e., the dynamic power splitting 
receiver was studied in [12]. The extension of [10] to the 
scenario with imperfect channel state information (CSI) at 
the transmitter was studied in [13]. For multiple-input single¬ 
output (MISO) channels, the optimal beamforming designs 
for SWIFT systems with/without secrecy constraint have been 
investigated in [14,15], and the optimal transmission strat¬ 
egy maximizing the system throughput of MISO interference 
channel has been studied in [16]. Moreover, the application of 
RF energy transfer technique in cognitive radio networks with 
multiple antennas at the secondary transmitter was considered 
in [17]. Finally, cellular networks with RF energy transfer were 
considered in [8,18]. It is worth noting that all these prior 
works focus on the point-to-point communication systems. 

The RF energy harvesting technique also finds important ap¬ 
plications in cooperative relaying networks, where an energy- 
constrained relay with limited battery reserves relies on some 
external charging mechanism to assist the transmission of 
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source information to the destination [2]. As such, a number 
of works have exploited the idea of achieving simultaneous 
information and energy transfer in cooperative relaying sys¬ 
tems [4,9,19-21]. Specifically, [20] studied the throughput 
performance of an amplify-and-forward (AF) relaying system 
for both time-switching and power-splitting protocols and [ 21 ] 
considered the power allocation strategies for decode-and- 
forward (DF) relaying system with multiple source-destination 
pairs. More recently, the performance of energy harvesting 
cooperative networks with randomly distributed users was 
studied in [4,9]. However, all these works are limited to 
the single antenna setup and all assume an interference free 
environment. 

Motivated by this, we consider a dual-hop AF relaying 
system where the source and destination are equipped with 
a single antenna while the relay is equipped with multiple 
antennas Q The energy constrained relay collects energy from 
ambient RF signals and uses the harvested energy to forward 
the information to the destination node. The power-splitting 
receiver architecture proposed in [10] is adopted. Specifically, 
we first study the performance of the multiple antenna relay 
system without CCI, which serves as a benchmark for the 
performance in the presence of CCI. Then, we present a 
detailed performance analysis for the system assuming a single 
dominant interferer at the relay. It is worth pointing out 
that, in the energy harvesting relaying system, while CCI 
provides additional energy, it corrupts the desired signal. In 
order to exploit CCI as a beneficial prospect, three differ¬ 
ent linear processing schemes, namely, 1) Maximum ratio 
combining/maximal ratio transmission (MRC/MRT), 2) Zero- 
forcing/MRT (ZF/MRT), 3) Minimum mean-square error/MRT 
(MMSE/MRT) are investigated. 

The main contributions of this paper are summarized as 
follows: 

• For the scenario without CCI, we derive an exact outage 
expression involving a single integral, and a tight closed- 
form outage probability lower bound. In addition, we 
present a simple high signal-to-noise ratio (SNR) ap¬ 
proximation, which reveals that the system achieves a 
diversity order of N, where N is the number of relay 
antennas. A new tight closed-form upper bound for the 
ergodic capacity is also derived. Finally, the optimal 
power splitting ratio minimizing the outage probability 
is characterized. 

• For the scenario with CCI, we present tight closed- 
form outage probability lower bounds and capacity upper 
bounds for all three schemes. In addition, we also char¬ 
acterize the high SNR outage behavior and show that 
both the MRC/MRT and MMSE/MRT schemes achieve 
a diversity order of N, while the ZE/MRT only achieves 
a diversity order of — 1. Moreover, the optimal power 
splitting ratio minimizing the outage probability is stud¬ 
ied. 

• The presented analytical expressions provide an efficient 

'This particular system setup is applicable in several practical scenarios 
where two nodes (e.g., machine-to-machine type low cost devices) exchange 
information with the assistance of an advanced terminal such as a cellulai* 
base-station/clusterhead sensor. [5,22,23] 


means to evaluate key system performance metrics, such 
as the outage probability and ergodic capacity, without 
resorting to time-consuming Monte Carlo simulations. 
Therefore, a fast assessment of the impact of various 
key system parameters such as the energy harvesting 
efficiency 77 , the number of antennas N, the source 
transmitting power pi and the interference power pi on 
the optimal power splitting ratio is enabled. 

• Our results demonstrate that the CCI could be potentially 
exploited to significantly improve the system’s perfor¬ 
mance. However, the actual performance gain due to 
CCI depends heavily on the choice of linear processing 
schemes. It is shown that the MMSE/MRT scheme is 
always capable of turning the CCI as a desired factor, 
and can achieve higher performance gain when the CCI 
is strong. On the other hand, CCI is not always beneficial 
when the MRC/MRT and ZE/MRT schemes are used. 
The performance degrades significantly in the strong 
CCI scenario if the MRC/MRT scheme is applied. In 
contrast, a weak interferer degrades the performance of 
the ZE/MRT scheme, which on the other hand achieves 
almost the same performance as the MMSE/MRT scheme 
in the presence of strong CCI. 

The remainder of the paper is organized as follows: Section 
II introduces the system model. Section III investigates of the 
performance of the system without CCI. Section IV addresses 
the scenario with CCI. Numerical results and discussions are 
provided in Section V. Einally, Section VI concludes the paper 
and summarizes the key findings. 

Notation: We use bold upper case letters to denote matrices, 
bold lower case letters to denote vectors and lower case letters 
to denote scalars. ||h||^ denotes the Erobenius norm; E{a;} 
stands for the expectation of the random variable x\ * denotes 
the conjugate operator, while T denotes the transpose operator 
and f denotes the conjugate transpose operator; CA/(0,1) 
denotes a scalar complex Gaussian distribution with zero mean 
and unit variance; r(a;) is the gamma function; ^{a,b;z) 
is the confluent hypergeometric function [24, Eq. (9.210.2)]; 
Ky (x) is the u-th order modified Bessel function of the second 
kind [24, Eq. (8.407.1)]; Ei(a:) is the exponential integral 
function [24, Eq. (8.211.1)]; F (a,x) is the upper incomplete 
gamma function [24, Eq. (8.350.2)]; 2 Fi{a,b;c; z) is the 
Gauss Hypergeometric function [24, Eq. (9.100)]; ip (x) is the 
Digamma function [24, Eq. (8.360.1)]; G™," (■) is the Meijer 
G-function [24, Eq. (9.301)] and (•) denotes the 

generalized Meijer G-function of two variables [25] which can 
be computed by the algorithm presented in [26, Table II]. 

H. System Model 

We consider a dual-hop multiple antenna AE energy har¬ 
vesting relaying system as shown in Eig. |l(a)| where both 
the source and the destination are equipped with a single 
antenna, while the relay is equipped with N antennas [5]. The 
source sends information to the destination through an energy 
constrained relay node. Throughout this paper, the following 
assumptions are adopted: 1) It is assumed that direct link 
between the source and the destination does not exist due 
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to obstacles and/or severe fading. 2) The channel remains 
constant over the block time T and varies independently and 
identically from one block to the other, and has a Rayleigh 
distributed magnitude. 3) As in [27-29], no CSI is assumed 
at the source, full CSI is assumed at the relay, and local CSI 
is assumed at the destination. 



(b) 


Fig. 1: (a) System model: S, R and D denote the source, 
relay and destination node, respectively, (b) Block diagram 
of the relay receiver with the power splitting protocol. 


We focus on the power splitting protocol proposed in [10]. 
Specifically, the entire communication consists of two time 
slots with duration of ^ each. At the end of the first phase, 
each antenna at the relay node splits the received source signal 
into two streams, one for energy harvesting and the other for 
information processing as depicted in Fig. |l(b)| As in [12,21], 
we consider the pessimistic case where power splitting only 
reduces the signal power, but not the noise power. Hence, our 
results provide a lower bound on the performance for practical 
systems. We now consider two separate cases depending on 
whether the relay is subject to CCI or not. 


A. Noise-limited Case 

Let 9 denote the power splitting raticJl, then the signal 
component at the input of the information receiver is given 
by 

yr = \J{^- 6»)Ps/d[hiX + tlr, (1) 


where Pg denotes the source power, hi is an W x 1 vec¬ 
tor with entries following identically and independently dis¬ 
tributed (i.i.d.) CAf{0, 1), di denotes the distance between 
the source and the relay, r is the path loss exponent, x is 
the source message with unit power, is an W x 1 vector 
and denotes the additive white Gaussian noise (AWGN) with 

Ejn^nJ;} = WqI. 

At the end of the first phase, the overall energy harvested 
during half of the block time can be expressed as 


Eh = 


r,ePg 


|hi| 


T 


dl """''^’ 2 ’ 
where p denotes the RF-to-DC conversion efficiency. 


( 2 ) 


B. Interference plus Noise Case 

We assume that the relay is subjected to a single dominant 
interferer and AWGN while the destination is still corrupted by 
the AWGN onlyH It is worth pointing out the single dominant 
interferer assumption has been widely adopt in the literature, 
see [31,32] and references therein. Moreover, such a system 
model enables us to gain key insights on the joint effect of 
CCI and multiple antennas in an energy harvesting relaying 
system. 

In such case, the signal at the input of the information 
receiver at the relay is given by 

= Y^(l - e)Ps/dl\iix + \J{1- 9)Pj/d}laiSi + n,., 

( 3 ) 

where Pj is the interference power, di denotes the distance 
between the interferer and the relay, si is the interference 
symbol with unit power, and h/ is an x 1 vector with 
entries following i.i.d. 1). 

Also, according to [11], at the end of the first phase, the 
overall energy harvested during half of the block time E is 
given by 

F;, = p0(^||hi||^ + ^||h,||^)|. (4) 

For both cases, during the second phase, the relay transmits 
a transformed version of the received signal to the destination 
using the harvested power. Hence, the signal at the destination 
can be expressed as 

Vd = \J^/dlh2Wyr -b Ud, (5) 

where h 2 is a 1 x W vector and denotes the relay-destination 
channel and its entries follow i.i.d. CAf{0, 1), d 2 denotes the 
distance between the relay and the destination, Ud is the 
AWGN at the destination with E{n'^nd} = Nq, W is the 
transformation matrix applied at the information receiver at 
the relay with E{||Wyj.||^} = Obviously, the performance 
of the system depends on the choice of W, which will be 
elaborated in the ensuing sections. 


^The optimality of uniform 9 can be established by using similar methods 
as in [12], 


^The scenaiio where the relay and the destination experience different 
interference patterns will occur in frequency-division relaying systems [30]. 
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III. The NOISE-LIMITED Scenario 


which can be evaluated as 


In this section, we consider the scenario where the relay 
is corrupted by AWGN only. In such case, it can be shown 
that the optimal transformation matrix W has the following 
structure; 


W = a; 




12 | 1 f 


|hi| 


( 6 ) 


where u! is the power constraint factor, i.e., the information 
receiver first applies the MRC principle to combine all the 
signals from N antennas, and then forward the signal to the 
destination by using the MRT principle. To guarantee the 
transmit power constraint at the relay, w can be computed 
as 


(jj 


2 


{i-e)Ps 

di 


_ 

l|hi||^ + Nq 


(7) 


where is the available relay power. Since the relay commu¬ 
nicates with the destination for half of the block time we 
— II hi 11^ . Hence, the end-to-end SNR 

of the system can be expressed as 


have = ^ = 2 ^ ||hi| 


UJ 


Ic 


hi 


7 = 


,2 (i-e)p. 

Ip drdi 


|h 2 ||" ^ + 


7)e(i-6>)p^ 


|h 2 


Ip I 


P\\f 


ff^l|h2f,.||hi 


Ip- 


(i-g)pi 

dT 


+ 1 


( 8 ) 


where pi is defined as pi = Ps/Nq. 

In the following, we give a detailed performance analysis 
in terms of the outage probability and ergodic capacity. In 
addition, the optimal 9 minimizing the outage probability is 
investigated. 


A. Outage Probability 

The outage probability is an important performance metric, 
which is defined as the instantaneous SNR falls below a pre¬ 
defined threshold 7 th. Mathematically, outage probability can 
be expressed as 

f’out = Prob(7 < 7th). (9) 


Theorem 1: The outage probability of the multiple antenna 
energy harvesting relaying system can be expressed as 


Tout — 1 ~ 




r,N-l 


r(iv) r(iv) 


^dx. 


( 10 ) 


P^PlTth 

dTdX ■ 


vhnrn n — r _ „ _ ’)S(l-«')p? j _ 

Where a - --, b - 7 th, c - , d- 

Proof: Substituting ([ 8 ]l into (|9]l, the outage probability of 

the system can be expressed as 


Pout= (11) 

Prob(||h 2 ||^ (c||hi||^ - d||hi||p) < (a ||hi||p-f 6 ^) , 



Since the squared Frobenius norm of a complex Gaussian 
vector is Chi-square distributed, ||hi||^ and ||h 2 ||^ are i.i.d. 
gamma random variables. After some simple algebraic manip¬ 
ulations (doll is obtained. □ 

Theorem[T]presents the exact outage probability of the of the 
system with arbitrary number of antennas. For the special case 
with a single antenna at the relay. Theorem [T] reduces to the 
result derived in [20, Proposition 3]. To the best of the authors’ 
knowledge, the integral in (O does not admit a closed- 
form expression. However, this single integral expression can 
be efficiently evaluated numerically using software such as 
Matlab or MATHEMATICA. Alternatively, we can use the 
following closed-form lower bound for the outage probability, 
which will be shown to be tight across the entire SNR range 
in the Section |V] 

Corollary 1: The outage probability of the multiple antenna 
energy harvesting relaying system can be lower bounded as 


plow _ -| _ 

-^OUt ~ 




N-1 N-1 

It. 

i—Q j—Q 




N-1 

j 








(13) 


Proof: See Appendix |T] □ 

While Theorem [T| and Corollary [T] are useful to study the 
system’s outage probability, the expressions are in general too 
complex to gain insight. Motivated by this, we now look into 
the high SNR regime, and derive a simple approximation for 
the outage probability, which enables the characterization of 
the achievable diversity order. 

Theorem 2: In the high SNR regime, i.e., pi oo, the 
outage probability of the multiple antenna energy harvesting 
relaying system can be approximated as 


poo 

-^OUt 


( 77 ) / 1 

r(iv + i) i (i_p) 


ln((l-0)pi)-ln(d[7th)-C/d: 


r(iv) 


N\ 


(14) 


where C is the Euler-Mascheroni constant [24, Eq. (9.73)]. 

Proof: See Appendix HIl □ 

We observe that the system achieves a diversity order of 
N, which is the same as the conventional case with constant 
power relay node [22]. However, we notice that Pout decays 
as p]”^lnpi rather than pf^ as in the conventional case 
[22]. This important observation implies that, in the energy 
harvesting case, the slope of Pout converges much slower 
compared with that in the constant power case. Please note 
that similar observations have been made in prior work [ 21 ]. 
The possible reason is that, in SWIPT systems, the available 






















5 


transmit power at the relay is a random variable, which results 
in higher outage probability compared to the conventional 
constant relay power case. 


B. Ergodic Capacity 

Noticing that the end-to-end SNR given in ([ 8 ]l can be 
alternatively expressed as 


7i72 

7i + 72 + 1 ’ 


(15) 


where 71 = ||hi||^ and 72 = ^ ||h 2 ||^ II hi II f- The 

ergodic capacity is given by 





7172 \ 

7i + 72 + 1 / 


(16) 


Unfortunately, an exact evaluation of the ergodic capacity is 
generally intractable, since the cumulative distribution function 
(c.d.f.) of dSll can not be given in closed-form. Motivated by 
this, we hereafter seek to deduce a tight bound for the ergodic 
capacity. 

Starting from (fTbl l. the ergodic capacity can be alternatively 
expressed as 




(1 -b7i) (1 -(-72) \ 

1 + 71+72 ) 


(17) 


where [log 2 (1 + 7 ^)], for i S {1,2}, and 

= 5 E [log 2 (1 + 71 + 72 )]. Using the fact that / (x, y) = 
log 2 (1 + + e^) is a convex function with respect to x and 

y, we have 

C-y^ > ^log2 (1 + . (18) 

We now establish the ergodic capacity upper bound of the 
system using the following theorem: 

Theorem 3: The ergodic capacity of the multiple antenna 
energy harvesting relaying system is upper bounded by 


61,0 = 


AT 1 


P 21n2 
1 


A;=0 


di 


N-1 




21n2r(N) ^ m! 

^ ' m—O 


(1-0) Pi 

didi 


r -k 


di 


ri0pi 


3,1 


Gt: 


dldl 


^log 2 


(1 - d)Pl MN) 

dl 


rjBpi 

yOpi 

dldl 


’(1-0) Pi 


— m 

— m, AT—m,0 




(19) 

□ 


Proof: See Appendix Hill 
Theorem [3] presents a new upper bound for the ergodic 
capacity of the system, which is quite tight across the entire 
SNR range as shown in the section [V] hence, providing 
an efficient means to evaluate the ergodic capacity without 
resorting to Monte Carlo simulations. In addition, as we show 
in the next subsection, it enables the study of the optimal 
power splitting ratio. 


more transmission power at the relay, which benefits the relay- 
destination transmission. Nevertheless, a large 0 also dete¬ 
riorates the quality of the source-relay transmission. Hence, 
there exists a delicate balance, which we now investigate. For 
tractability, we only focus on the outage performance in the 
high SNR region, and the impact of 0 on the ergodic capacity 
will be numerically illustrated in the Section |V] 

Starting from the high SNR approximation of Pout in (O, 
the optimal 0 , which is the solution of the optimization 

problem min Pout, can be obtained by solving the equivalent 
0<6I<_L 

problem in (EOb shown on the top of the next page. 

Proposition 1: The optimal 0 is the root of the following 
polynomial 


ai0jv+i _ _ 0)^+1 _ _ 0)^ 

-di(l-0)^+Mn(l-0) =0, (21) 


where ai = N, by = 


Cl = 


? 7 «f(Ar)’ 


dl = 


Afd, 


and 0 < 0 < 1 . 


Proof: It is easily to prove that, when pi — 00 , there is 
only one root (denoted by 0 *) on the interval of ( 0 , 1 ) for the 
equation f'{9) = 0, and we can also note that /'(O) = —00 
and /'(I) = + 00 . Due to the continuity of we have 

fie) < 0 , 0 € ( 0 , 0 *) and /'( 0 ) > 0 , 0 € ( 0 *, 1 ), which 
means that /( 0 ) first decreases as 0 from 0 to 0 * and then 
increases as 0 from 0* to 1. Therefore, the global minimum 
of /( 0 ) can be obtained by solving f ( 0 ) = 0 . □ 

Due to the presence of logarithmic function In (1 — 0), a 
closed-form expression for the root of (I 2 TI 1 can not be obtained. 
However, it can be efficiently solved numerically. 


TV. The Interference plus Noise Scenario 

We now assume that the relay is subject to the influence 
of a single dominant interferer. In the presence of CCI, the 
optimal relay processing matrix W maximizing the end-to-end 
signal-to-interference-and-noise ratio (SINR) of the system is 
the solution of the following optimization problem: 

^^|h2Whir 

ITLclX ^ _____ 

s.t. E{||Wy,||^} = Pr = ye llhill^ + ^ llhill^^ . 

( 22 ) 

Due to the non-convex nature of the problem, a closed-form 
solution for W is hard to And. Hence, in the following, we 
consider three heuristic two-stage relay processing strategies 
proposed in [23], i.e., the matrix W admits the rank-1 structure 

ht 

W = Wi, where wi is a 1 x TV linear combining vec¬ 

tor, which depends on the linear combining scheme employed 
at the relay and will be specified in the following subsection. 


C. Optimization of the Parameter 6 in High SNR Value 

The right selection of the power splitting ratio 0 is crucial 
for the system’s performance. A high value of 0 could provide 


A. MRC/MRT Scheme 

For the MRC/MRT scheme, Wi is set to match the Arst hop 
channel given in (| 6 l). To meet the transmit power constraint at 
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min / iO) 
0<9<1 ^ 


1 ln((l - 9) pi) - lnd[7th " C / 

(i-< r(iV) 


( 20 ) 


the relay, the power constraint factor should be given by 

Pr 

-, (23) 


a;2 = 


(i-e)Ps 

di 


I hi 


,2 (l-e)Pj Ihjhj 


iihiiii: 


Nn 


where Pr = = 


. Therefore, 


»7^'(^l|hi||^ + ^||h/||p^ 

the end-to-end SINR of the MRC/MRT scheme can be ex¬ 
pressed as 


,^IVIRC,^MRC 

MRC _ 7/1 7/2 _ 

2'^ 7/T"+7/T" + 1’ 


(24) 


where 


MRC 

7/1 


{i-e)pi I 

—31-1 


^1 IIf 


(i-e)p, 

—3T- 


iij hj 


MRC 

7/2 


--I-1 


and Pi is defined as 


7) Outage Probability: Since the exact analysis appears to 
be difficult, in the following we focus on deriving an outage 
lower bound and a simple high SNR outage approximation. 
According to [22,23], the end-to-end SINR in (l24li can be 
tightly upper bounded by 


MRC ^ up 

7/ < 7/ = mm 


(7. 


MRC 


MRC 

7/2 


(25) 


the outage probability of the MRC/MRT scheme is lower 
bounded by 

P|„^f^ = Prob(77<7th). (26) 


We have the following key result: 

Theorem 4: If pi ^ p/0 the outage probability of the 
MRC/MRT scheme can be lower bounded as 

r>LMRC -I ttiMRCttiMRC /o'TX 

^/out —-L — ^2 1 

with 

pMRC ^ V 

^ (m-n)! Vdjpi-f d[p/7th/ 

and can be expressed as in (|28] | shown on the top of 

the next page. 

Proof: See Appendix IIVI □ 

While Theorem|4]is useful for the evaluation of the system’s 
outage probability, the expression is too complex to yield 
much useful insights. Motivated by this, we now look into 
the high SNR region, and derive a simple approximation for 
the outage probability, which enables the characterization of 
the achievable diversity order of the system. 


^For mathematical tractability, we only provide the result for the general 
case where the signal from the source and the CCI have different power, i.e., 
pi ^ p/. But the result for the special case pi = pi is much more easier 
and can be obtained in a similar way. 


Theorem 5: In the high SNR region, i.e., pi 00, the out¬ 
age probability of the MRC/MRT scheme can be approximated 

aF 


dMRC 

/out 





I- 9 J ^dT{N-n)\ 

/ n=0 ^ ^ ’ 


N-1 
2=0 


2-Fi N,2N-i-l-,2N-i-,l- 


\ 

^?Pi J 


2N-i-l 


{r]e)^r{N + 1)T{N) 


Proof: See Appendix [V] 


/ 

(29) 

□ 


For the special case where the relay is equipped with a single 
antenna, i.e., TV = 1, with the help of [24, Eq. (9.121.6)], ( i29b 
reduces to 


pMRC 

^lout 



EL 

d} 


^2 (In In If: 


r]6 


dljth 

Pi 


(30) 


Theorem |5] indicates that a full diversity order of N is 
still achievable in the presence of CCI for the MRC/MRT 
scheme. Moreover, from (i30l l, we see that the effect of CCI 
could be either beneficial or detrimental, depending on the 
relationship between p/, df, dj V and 6, i.e., when — 

positive, the CCI is detrimental, while when 
d.j) negative, the CCI becomes beneficial, 
witich suggests that, in wireless powered relaying systems, 
CCI could be potentially exploited to improve the perfor¬ 
mance. 


2) Ergodic Capacity: Utilizing similar techniques as in 
Section IIII-BI we establish the following ergodic capacity 
upper bound: 

Theorem 6: If pi ^ p/, the ergodic capacity of the 
MRC/MRT scheme is upper bounded by 


/^MKC I 

O/Lip = O^MRC + O^MRC — 

+ + (31) 

where C^mrc, C^urc, E (ln7^j^'^‘') and E are given 

by (I32I 1 - (|35]l shown on the next page. 

Proof: See Appendix IVII □ 

3) Optimal 0 Analysis: We now study the optimal value of 
9 minimizing the outage probability. Based on the high SNR 
approximation for in ( i29] l, the optimal 9 can be found 

as: 

Proposition 2: The optimal 0 is a root of the following 


^It is worth pointing out that the result in Theorem [^ holds for all cases 
whether the signal power and the CCI power is equal or not. 
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^MRC _ 2^1 


iNtjNt ^ 


PlP^ 


•Y ■ ^ ii-N- j) 


■2. 


(A^-s)!(s-l)! \pi Pi 


l — N—sN—l / IT \ 

J_ / »2 7th A 

^ ml\p9 J 


N+l-s 


m—0 



2<-<- ^ UUl (l-N- j) 
PiP? h {N-s)\{s-l)\ 


dld}jth \ ' 

pOpi ) 


^m+s — N—l I 


I dldJjth 

pOpi 


(28) 


C^MRC 


(1 - 6) pi y^ y^ ((1 - 0) piY" 1 , 14 , 1,1 

2dTln2 ^ dli'^im — n)\n\ 

^ 771=0 71=0 


(1-0)P1 

(i-e)p, 


771+1 
0 ; —71 

0;0 


(32) 




,MRC TtT TV T 

pfpf21n2 


AT 

E 

S = 1 


{N - s)! (s - 1)! \pi 


G+ 


fdldl 
3 Vl^Pi 


S-7V-1 

S —A/” —1,771 + S —iV—1,0 


pi 


I-N-S N-1 


1 

to ! 


Af+l-s 


pf+21n2 


771 = 0 

TTS-1 

E 


S=1 

N-1 


^ - d) (dl d] 

{N - s)! (s - 1)1 El pi 


771=0 


1 

to ! 


N+l-s 


Gt 


( did} 

\p0pi 


1-N-s 


s-N-1 

s —N — 1,m-\-s — N—1,0 


(33) 


E (In7/T^) = ln((l - 0) pi) - lnd[ + ^ (1) - eT^ 

N—1 m 

EE 


d} 


1,1 

0,0,1 


0 - 0 ) PI 
{{i-0)pN-"‘ 

(tI —77l)' 


=i„=o {m - n)\d\ 


tE (to) il' 


f TO, 


m — n; 


0 - 0 ) PI 


(34) 


polynomial 


N-1 


E-4(n)(l- + 


i-N-l 


B 


71 = 0 


pw+i 


pj 


= 0 , 


B 


where A{n) = (n-u-i)v 

,Nt n-1 2Fi( N,2N-i-l-2N-i-l-^I^] 

d2 Y- (N-l\( ..N-i-l ^ A _ 

i7«r2(Ar) + ( z )V I) 

and 0 < 0 < 1. 


2N-i-l 


In the special case of = 1, the optimal solution can be 
given in closed-form as follows: 


(36) 


^opt _ 
■^MRC ~ 


pi (in pj —In pi —In -t-ln d\ ) 


ri[dlpi-d^jPi) 


dj pi (in pj- —In pi —In -t-ln d \) 

p(dlP/-<iJpi) 


1-f 


(37) 


This simple expression is quite informative, and it can be 
observed that the optimal 9 in (l37l i is a decreasing function 
of p and Pi, and an increasing function of pi, which can be 
explained as follows: 


Proof: The result is derived by following the same steps 
as in the proof of Proposition [T] □ 


As p increases, more transmission power can be collected 
at the relay, hence the bottleneck of the system perfor- 


jWt N -prs-l Q ^ _ 1 

E (In7/T^) =\iip9-\ndl+^{N) + ^Y2- 

Pi s=i 


(5-1)! 

N T-r«-i 


PI Pi 


('(/' (iV — s -b 1) + In Pi — In c?i 


Pi 


(5-1)! 


VPI 



1-S 


(A^ — 5 -h 1) + In p/ — In d}). (35) 


















































mance lies in the SINR of the signal at the input of the 
information receiver. As a result, we should choose a 
smaller 0 to improve the first hop performance. 

• A large pi provides more energy, while at the same time 
reduces the SINR of the first hop transmission. Hence, a 
smaller 9 should be chosen to compensate the loss of the 
SINR. 

• For large pi, in general the first hop transmission quality 
is quite good, hence, it is beneficial to have more energy 
at the relay, i.e., a larger 0 is desirable. 


B. ZF/MRT Scheme 


For the ZF/MRT scheme, the relay utilizes the available 
multiple antennas to completely eliminate the CCI. To ensure 
this is possible, the number of the antennas equipped at the 
relay should be greater than the number of interferers. Hence, 
for the ZF/MRT scheme, it is assumed that TV > 1. According 
to [23], the optimal combining vector Wi is given by Wi = 

—i==, where P = Iat — h/fhth/) ht. Therefore, the 

yhJPhi V / 

end-to-end SINR of the ZF/MRT scheme can be expressed as 


iF = 




7?r + 7E + i’ 


(38) 


where 

f (tiii>.iiT 


fllh/ 


hlPhi 




(l-g)pl 

dJ ' 


7/2 


1) Outage Probability: We first present the following out¬ 
age lower bound: 

Theorem 7: If p\ ^ pi, the outage probability of the 
ZF/MRT scheme can be lower bounded as 


pLZF _ 1 pZFpZF 
^/out — ^2 I 


(39) 


Where ^ = F^^^. 

Proof: According to [23], the c.d.f. of is given by 

dT X 


F W-1, 


(^) = 1 “ 


(i-e)p 


r) 


r{N-i) 


(40) 


Then, the desired result can be obtained by following the 
similar lines as in the proof of Theorem |4] □ 

To gain further insights, we now look into the high SNR 
region, and present a simple and informative approximation 
for the outage probability. 

Theorem 8: In the high SNR region, i.e., pi —F oo, the out¬ 
age probability of the ZF/MRT scheme can be approximated 
as 


pZF 


rf[7th 

(7V-1)! V(l-0) Pi 


N-l 


(41) 


Proof: With the help of the asymptotic expansion of 
incomplete gamma function given in [24, Eq. (8.354.2)], it 
is easy to note that the c.d.f. of can be approximated as 


F zf 
<11 


(x) 


d^x 


(iV-1)! V(l-0)pi 


N-l 


(42) 


Then, utilizing (i42l i and following the similar lines as in the 
proof of Theorem |3 we can obtain 


F 


ZF 

/out ' 


rf[7th 

(TV-1)! V(l-0 )Pi 


N-l 


f d^ ^2 Tth \ 

V ) 


N 


r(TV + i)r(TV) f 


N-l 

E 


i=0 


A^- 1 




2^1 


fTV, 2TV - i - 1; 2TV - T; 1 - ^ 

V ’ ’ ’ d,iPi 

2N -i -I 


• (43) 


The desired result follows by noticing that the second term is 
negligible compared with the first term in (021. □ 

Theorem |2 indicates that the achievable diversity order of 
the ZF/MRT scheme is TV — 1. Compared with the MRC/MRT 
scheme, the ZF/MRT scheme incurs a diversity loss of one. 
This is an intuitive and satisfying result since one degree of 
freedom is used for the elimination of the CCI. 

2 ) Ergodic Capacity: We now look into the ergodic capac¬ 
ity of the system, and we can establish the following upper 
bound of the ergodic capacity: 

Theorem 9: If pi pi, the ergodic capacity of the ZF/MRT 
scheme is upper bounded by 

CFp = - ilog2 (l + ^ 

(44) 


where C'-.zf 

'll 


e {i-®)Pi 
2 In 2 


N-2 / 

s ( 


( 

A 


V(1 

-6)pi 

; ^ (i-e)PiJ 


E (In7/0 = In ((1- 6*) Pi) - lnd[ + ip {N - 1), C.^zf = 
C^uRc and E F^lF) = E (in7^2*^^). 

Proof: With the help of the c.d.f. of 7jf given in (i40l) 
and following the similar lines as in the proof of Theorem |6] 
yields the desired result. □ 

3) Optimal 9 Analysis: We now study the optimal 9 
minimizing the outage probability. Based on the high SNR 
approximation for 02 , the optimal 9 can be found 

as: 

Proposition 3: The optimal 0 is a root of the following 
polynomial 


Ai 

{1-9)N 


Bi 

QN+l 


= 0 , 


(45) 


where Ai = 

V'' i\N—i — l 

i7«r2(Ar)pi \ i 


1 K? _ 

(Ar-2)!’ “ 

2Fi(^JV.2Ar-*-l;2Ar-*;l-^) 
2Ar-i-l 


and 0 < 0 < 1. 

Proof: The result is derived by following the same steps 
as in the proof of Proposition [T] □ 


C. MMSE/MRT Scheme 

The ZF scheme completely eliminates the CCI at the 
relay, which however may cause an elevated noise level. In 
contrast, the MMSE scheme does not fully eliminate the CCI, 
instead, it provides the optimum trade-off between interference 
suppression and noise enhancement. According to [23], Wi is 
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given by 


Wi = h| [ hih| + h/h 


t 


,t 


(46) 


Therefore, the end-to-end SINR of the MMSE/MRT scheme 


can be expressed as 


,^MMSE 


„,MMSE„,MMSE 

In 112 

„,MMSE I ,^MMSE _ 
'll 112 


„,MMSE _ vS 

T/2 - T 


dT 


+ ^llh 

I 


/IIf 


(i-e)p/ 


(47) 


I and 




where R = h/h 

f (t iiCi 

1) Outage Probability: 

Theorem 10: If pi ^ pi, the outage probability of the 
MMSE/MRT scheme can be lower bounded as 


pLMMSE _ 1 7 /MMSEz?MMSE 
^/out — ^ 1 


where i^MMSE 


<^1 Tth 

e ( 1 -»)P 1 (l-e)pj 
d-T{N) 




<il7th ' 

PI, 


r(Ar) 




(48) 


and 


^MMSE _ ^MRC 

Proof: According to [23] we know that the c.d.f. of 

^MMSE jg gjygjj ]3y 


TLmmse 

'/I 


(a;) = 1 - 


("•ifSk) 


e (i-»)pi (1 — 9) Pi 


dfx 


(1-0) PI 


r(iv) 

N 

2 F 1 


<PjT {N) 


,l;2;-/l«x 


d}pi 


(49) 


Then, following the similar lines as in the proof of Theorem 
|4] we can obtain the desired result. □ 

To gain further insights, we now look into the high SNR 
region, and present a simple approximation for the outage 
probability. 

Theorem 11: In the high SNR region, i.e., pi —00, 
the outage probability of the MMSE/MRT scheme can be 
approximated as 


pMMSE 



r (iv-f 1) r (/V) 



(i-g)p/ A 1 

d-r(iV) ) (i_0)^ 


N-1 


E 


TV- 1 
i 


I \ N—i — \ 

(-1) X 


2Fy (/V,2/V-/-1;2/V-/;1 
2N -i-\ 



(50) 


Proof: After some simple manipulation the c.d.f. of 
^^MSE approximated as 


F MMSE (x) 

III ' ' 



(1 -9)pi \ / d[7th N 
d-r(iV) J \ il - 9) pj 


(51) 


Then, following the similar lines as in the proof of Theorem 
|5] we can obtain the desired result. □ 

Theorem [TT] indicates that the MMSE/MRT scheme 
achieves a diversity order of N, the same as the MRC/MRT 


scheme. 

A close observation of (l29l l. (l43T l and (fSOl l reveals that the 
difference among all three schemes only lies in their first 
terms, which can be expressed as follows: 


N 

omrc = ^ 

n—0 


d^'^(N-n)V 


= (iv^’ 

1 , {i-e)pi 


(52) 


It can be easily observed that ommse is strictly smaller 
than omrc, since ommse only includes the first two terms of 
omrc- As such, we conclude that the MMSE/MRT scheme 
always achieves a strictly better outage performance than 
the MRC/MRT scheme due to the higher array gain. For 
the ZF/MRT scheme, although a diversity loss leads to its 
inferior performance in the high SNR region, it should be 
noted that ozf is generally smaller than omrc, which means 
that the ZF/MRC scheme has a larger array gain than the 
MRC/MRT scheme. Therefore, in the low SNR region, the 
ZF/MRT scheme may achieve better outage performance than 
the MRC/MRT scheme. 

2) Ergodic Capacity: 

Theorem 12: If pi 7/ pi, the ergodic capacity of the 
MMSE/MRT scheme is upper bounded by 


Flup = C^MMSE 


■ C..„MMSE — 
>12 




where 




21n2 


(1 — 9) p^pi 1 , 1 , 2 , 1,1 


21n2 r(iV)df d[^l’[l:2],0,[l:2] I 


(i-e)pi 

dT 


Ar+2 ^ 

0;(-2,-l) 
0;(-i-2), 


(54) 


E (In7/T®^) =fiN)+ In ((1 - 0)pi) - lnd[- 


((1 - d)pif ^ 1^3 

dfr(/v) 


( a-d)pi 
V d} 


-N,-2, 

- 1,-2 



(55) 


and C.^MMSE = C^wRc as well as E (in 7^''^^^) = E (in 7}^'^^). 

Proof: With the help of the c.d.f. of gjygjj jjj 

and follows the similar lines as in the proof of Theorem |6] 
yields the desired result. □ 

3) Optimal 9 Analysis: We now study the optimal 9 
minimizing the outage probability. Based on the high SNR 
approximation for in (l50t . the optimal 9 can be found 

as: 

Proposition 4: The optimal 0 is a root of the following 
equation 


1 I {N-l)pi B 

(1 - 0)^+ir(Af) d}(i - 0)^r(w) 0^+1 

where B have been defined in (1361) and 0 < 0 < 1. 


(56) 
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Proof: The result is derived by following the same steps 
as in the proof of Proposition [T] □ 

V. Numerical Results and Discussion 

In this section, we present numerical results to validate the 
analytical expressions presented in Section IV, and investigate 
the impact of various key system parameters on the system’s 
performance. Unless otherwise specified, we set 7th = 0 dB, 
rj = 0.8, 6 = 0.5, Pi = 9.5 dB, t = 2 and di = d 2 = dj = 1. 

A. Effect of Multiple Antennas 



(a) Outage probability 



(b) Ergodic capacity 

Fig. 2: Impact of N on the system performance. 

Fig. |2] illustrates the impact of antenna number N on the 
outage probability and ergodic capacity. It can be readily 
observed from Fig. |2(a)| that for all the three considered 
schemes, the proposed lower bounds in (l27l i. ( |39] | and (l48l l 
are sufficiently tight across the entire SNR range of interest, 
especially when N is large, and become almost exact in the 
high SNR region, while the high SNR approximations in (|2^ . 
(EB and ([soi l work quite well even at moderate SNR values 
(i.e.. Pi = 20 dB). In addition, we see that both the MRC/MRT 
and MMSE/MRT schemes achieve the full diversity order of 
N, while the ZF/MRT scheme only achieves a diversity order 
of TV — 1, which is consistent with our analytical results. 
Moreover, the MMSE/MRT scheme always attains the best 
outage performance among all three proposed schemes, and 


the ZE/MRT scheme outperforms the MRC/MRT scheme in 
the low SNR region, while the opposite holds in the high SNR 
region. 

Erom Eig. |2(b)| we see that, for all three schemes, the 
proposed ergodic capacity upper bounds in (IbTI i. (l4^ and (l5^ 
are sufficiently tight across the entire SNR range of interest. 
In addition, we observe the intuitive result that increasing N 
results in an improvement of the ergodic capacity. Moreover, 
the MMSE/MRT scheme always has the best performance, 
while the ZE/MRT scheme is slightly inferior, and the perfor¬ 
mance gap between them disappears as N increases. On the 
other hand, the MRC/MRT scheme always yields the lowest 
ergodic capacity, and as N increases, the performance gap 
becomes more pronounced. 

B. Effect of CCI 



(a) Outage probability 



(b) Ergodic capacity 

Eig. 3: Impact of CCI on the system performance. 

Eig. [3 investigates the impact CCI on the system perfor¬ 
mance. The scenario without CCI is also plotted for compari¬ 
son. It can be readily observed from Eig. |3(a)| that the outage 
probability of the MRC/MRT scheme decreases slightly for 
smaller pj (i.e., pi < 0 dB), and then increases as the inter¬ 
ference becomes stronger. This phenomenon clearly indicates 
that the CCI can cause either beneficial or harmful effect on the 
system’s performance. This is because that CCI provides addi¬ 
tional energy but at the same time corrupts the desired signal. 
Eor the MRC/MRT scheme, when the CCI is too strong, the 
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disadvantage of the CCI becomes the dominant performance 
limiting factor. However, with sophisticated interference miti¬ 
gation schemes, e.g., the ZF/MRT and MMSE/MRT schemes, 
such undesirable effect could be eliminated. As shown in these 
two schemes, the outage probability decreases monotonically 
as Pi increases. Moreover, for the MMSE/MRT scheme, CCI 
is always desirable, while for the ZE/MRT scheme, whether 
CCI is beneficial or not depends on its power. 


C. Effect of the Distance 



(a) Impact of relay location, cZ 2 = 5 — di, dj = 3. 



(b) Impact of interferer location, d\ = 2, ^2 = 3. 


Eig. 4; Impact of distance on the system performance. 

Eig.lHshows the effect of the node distances on the system’s 
outage probability. Unlike the conventional dual-hop system, 
where it is in general desirable to place the relay in the middle 
of the source and the destination, Eig. |4(a)| indicates that in 
the energy harvesting scenario, the optimal relay location tends 
to be close to the source. This observation implies that, the 
quality of the first hop channel is more important than that of 
the second hop channel. This is quite intuitive since the quality 
of the first hop channel not only affects the received signal 
power at the relay but also determines the available power 
for the second hop transmission. As shown in Eig. |4(b)[ as 
the distance of the CCI increases, the outage performance of 
ZE/MRT and MMSE/MRT schemes deteriorates. In contrast, 
the outage performance of the MRC/MRT scheme improves. 
This is also intuitive, since increasing the distance reduces 


the received power at the relay, which in turn deteriorates 
the performance of the ZE/MRT and MMSE/MRT schemes, 
since strong CCI is desirable for the both the ZE/MRT and 
MMSE/MRT schemes as illustrated in Eig. [3] 

D. Effect of Power Splitting Ratio 9 



e 


(a) Outage probability 



(b) Ergodic capacity 

Eig. 5: The optimal power splitting ratio 9. 

Eig.|5]investigates the the impact of the power splitting ratio 
9 on the outage performance. We observe that there exists 
a unique 9 which gives the best outage or ergodic capacity 
performance. Eor all three schemes, we see a similar trends 
on the impact of 9, i.e., when 9 increases from zero to the 
optimal value, the performance improves; when 9 exceeds the 
optimal values, the performance deteriorates gradually. This 
phenomenon is rather intuitive, since performance of dual- 
hop systems is limited by the weakest hop quality. Moreover, 
we see that the optimal 9 is in general different for different 
schemes and performance metrics, as shown in Eig. |5(a)[ the 
MMSE/MRT scheme requires a smaller 9 compared with the 
MRC/MRT scheme, and the capacity optimal 9 is larger than 
the outage optimal 9 for the MMSE/MRT scheme. 

E. Effect of Key System Parameters on the Optimal 9 

Eig. |6] examines the effect of various key system parameters 
such as p, N, Pi and pj on the choice of optimal 9. Specifi¬ 
cally, Eig. |6(a)| illustrates the effect of p, and we can see that. 
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the outage optimal 0 is a decreasing function of ij. A large 
r] implies higher energy conversion efficiency, which in turn 
suggests that less portion of the signal is needed for energy 
harvesting, hence, a smaller 0 is required. A similar trend is 
observed in Fig. |6(b)| on the impact of N. As N increases, the 
additional antennas improve the energy harvesting capability, 
i.e., more energy could be harvested, which implies that the 
optimal 9 should decrease. The effect of pi is shown in Fig. 
|6(c)| it is interesting to see that, for the MRC/MRT and 
MMSE/MRT schemes, the optimal 6 is an increasing function 
of pi, while for the ZF/MRT scheme, the optimal 9 increases 
first along pi and decreases when pi exceeds certain value. 
Finally, Fig. |6(d)| investigates the effect of pi. For all three 
schemes, the optimal 0 is a decreasing function of pj. This 
is intuitive since the CCI serves as the energy source, when 
the CCI power increases, a smaller 0 is sufficient to fulfill 
the energy requirement at the relay. Similar trends could be 
observed for the capacity optimal 0, which are not presented 
here due to space limitation. 



(a) AT = 2, Pi = 10 dB, pj = 
0 dB 




(b) p = 0.8, Pi = 10 dB, Pi 
0 dB 



(d) N = 2, rj = 0.8, Pi 
15 dB 


Fig. 6: The impact of (a) rj, (b) N, (c) pi, (d) pi on the 
outage optimal 0. 


VI. Conclusion 

In this paper, we studied the performance of a dual-hop 
AF energy harvesting system with multiple antennas and 
CCI. Analytical expressions for the outage probability, ergodic 
capacity, as well as the diversity order were presented, which 
provide efficient means for the evaluation of the system’s 
performance. In addition, the optimal power splitting ratio 
minimizing the outage probability was analytically charac¬ 
terized while the capacity optimal power splitting ratio was 
studied numerically. Moreover, the impact of various key 
system parameters, such as p, N, pi and pi on the optimal 0 
were examined, which provided useful design insights on the 
choice of a proper power splitting ratio under different system 
configurations. 

Our results demonstrate that both the MRC/MRT and 
MMSE/MRT schemes achieve a full diversity of N while 


the ZE/MRT scheme only attains a diversity order of — 1. 
We showed that the CCI could be potentially exploited to 
significantly improve the system’s performance. With the 
MMSE/MRT scheme, the CCI is always a desirable factor, 
and the stronger the CCI, the better the performance. Never¬ 
theless, this is not the case for the MRC/MRT and ZE/MRT 
schemes, where the CCI could be detrimental. Eor instance, 
strong interference degrades the system performance of the 
MRC/MRT scheme, whereas the performance worsens in the 
presence of weak interference with the ZE/MRC scheme. 


Appendix I 

Proof of Corollary[T] 


We first notice that the end-to-end SNR of the system can 
be tightly upper bounded by 


7 < 




12|If 


I hi 




yOpi 

d^dl 


IhlllU^i^llhil 


(57) 


Hence, we get the following the outage probability lower 
bound; 


P'Z = Prob (llh^ll^ (c llhill^ - d llhill^) < a llhill^) , 

(58) 


which can be computed as 



Noticing that ||hi||^ and ||h2||^ are i.i.d. gamma random 
variables, we have 



r{N) T{Nf ^ 


(60) 


Then, making a change of variable cx — d = t, dbOl) can be 
alternatively written as 


plow _ -I _ 

-^OUt ~ ^ ~ 



r{N,a/t) 

TiN) 


{t + d)^ 

( 61 ) 


Invoking the series expansion of incomplete gamma function 

[24, Eq. (8.352.4)] and applying the binomial expansion 

, N-l 

(t + d) = Yl, d^~^~^, (1^ can be further 

i=o ^ 

expressed as 


plow _ -| _ 

-^OUt ~ ^ 


N 


o-d/c 


N-1 , N-l 

sis 

j—0 
roo 


N-l 




-.7-1 


X [ (62) 

Jo 


To this end, with the help of [24, Eq. (8.432.7)], the desired 
result can be obtained. 
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Appendix II 
Proof of Theorem[2] 

Starting from (|8ll, we observe that, as pi — oo, the end-to- 
end SNR can be tightly bounded by 


^ UD 

7 < 7 = min 


(l-g)Pi 

dl 


dIIf > 


■nOpi 

didi 


12|If llhlllj. 

(63) 


We now study the c.d.f. of 7 ^p. Noticing that 7*^^ = 

, we first 


tllhil 


p Y, where Y = min ((1 - 6») , ^ ||h2|l^ 


look at the c.d.f. of Y, which can be expressed as 
Fy (y) = Prob (^Wh^Wl < ^ , ||h2||^ < 


Prob [1-9 <y , ||h2||^ > 


r]6 


(64) 


with 


Pi = 



P2 


. ||2 ^ (l-e)dj 


Vi ii2 ^ 


9) ’ 


, y>l-9 
y <1-9 


P2 = 


Prob(^||h2|l^ > 
0 , 


(i-g)rf; 

r,e 


, y>l-9 
y <1-9 

Therefore, the c.d.f. of Y can be finally expressed as 

1 , ^ y>l-0 

12|If < ^) ,y<i-9 

y>l-9 
y <1-9 


(65) 


Fy (y) = 



( 66 ) 


Having obtained the c.d.f. of Y, we are ready to compute 
the c.d.f. of 7‘^P as follows; 

r°° f zd7\ 

F^up {z)= Fy [ — ] /iihiii^ (x) dx, (67) 
Jo \PlX/ 

which can be expressed as 


Fyup {z) = l- 


/ OO 


(N, 

y ’ r)6pix J 


Ig a 


r (N) P (N) 


dx. (68) 


(i-»)pi 

Now, applying the asymptotic expansion of incomplete gamma 
function [24, Eq. (8.354.2)] to (l68l l yields 


Fyup (Z) Ri 1 — 


(i--e)p 


1 fzdldl 
N\ \ri9pix 






TiN) 


-dx. 


(69) 


Please note, due to the omission of the higher order items of 
the asymptotic expansion of incomplete gamma function, the 
expression given in (l69l l is no longer a bound, but a very tight 
asymptotic approximation, and matches well with the exact 
value in the high SNR region, i.e., —?> Prob (7 < 7th). 

Pi —>-oo 

To this end, with the help of [24, Eq. (8.350.2)] and [24, Eq. 
(8.211.1)], we obtain the following closed-form expression for 


poo 

-^OUt 

T{N) m{N-iy. V v9pi J ■ 

(70) 

Einally, applying [24, Eq. (8.214.1)] and [24, Eq. (8.354.2)] 
yields the desired result. 


Appendix III 
Proof of Theorem[3 


The ergodic capacity can be upper bounded by 


Cup = 


■ Cry^ 


1 


log2 ( 


1 + e 


,E(ln7i) I E(ln72) 


). (71) 


Note that is the ergodic capacity of the SIMO Rayleigh 
channel, which has been given in [37] as 


~ 21n2 


N-1 

E 


dl 


(1-0) Pi 


;(i-e)pi r I —k 


dl 


’(1-0) Pi 

(72) 


and is the ergodic capacity of the SIMO keyhole channel, 
which has been given in [38] as 


a 


72 


1 1 

21n2r(Af) 


N-l 


E 

m—O 



dldl 


p0pi 


— m 

— m,Ar—m,0 


(73) 


Next, we observe that E(ln7i) = In + 

E An||hi||^^ and E(ln72) = In (7^) + E (^ln||hi7^ -f 
E (in ||h2||^ ). It is easy to show that 

E(ln||hi||^) =E(ln||h2||^) = 9/9 ( N ). (74) 

To this end, pulling everything together yeilds the desired 
result. 


Appendix IV 
Proof of Theorem|4] 

Prom ( i26l l. the outage lower bound can be evaluated as 

Pi^ut^ = Prob (7rT^ < 7th) + Prob (7r2^^ < 7th) 

- Prob (7n'^^ < 7th and 7^2'^^ < 7th) 

r: Prob (7^1'^^ < 7th) + Prob (7^2'^^ < 7th) 

- Prob (77T^ < 7th) Prob (7/T^ < 7th) • (75) 

In general, 7^1'^*" and 7^'^*" are not independent. 
However, through Monte Carlo simulations, we 
observe that as long as p/ is close to pi, the term 
Prob (777'^'” < 7th and 7^^'^*" < 7th) can be safely 
approximated by Prob (7j^i'^^ < 7th) Prob (7)^'^'" < 7th) 
in the whole SNR region as shown in Pig |7(a)| As a 
matter of fact, the same approximation has been adopted 
in [39]. Therefore, the remaining task is to compute 
Prob (77^!'^^ < 7th) and Prob (77^2'^'= < 7th). 
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(a) Probability versus pi in dB. 



(b) Probability versus pi in dB. 


Fig. 7: Justification of the approximations employed in the 
proof of Theorem 4 and 5. 


(8.312.2)], we arrive at 


"^th AT 1 

F^urc (7th) = 1 - 


m 


a-o)pi 

1 


■E 

m—O 


d[7th 

(1-0) Pi 


X E 

^ (m - n)! yd}pi + d[p/7th 


(1 - 0) pipi 


n+1 


• (78) 


Similarly, the c.d.f. of 7^’^'” can be expressed as 
F.^mrc (7th) = Prob (7^2*^^ < 7th) = Prob |^||h2||E 

= (79) 


where Z = ■^ l|hi||^ + ^ ||h/||^, which is a sum of two 
independent gamma random variables. According to [35], if 
Pi Pi the probability density function (p.d.f.) of Z can be 
given by 


fz{x) = 


Pi P^ 


El 

PI 

TS-l 


Pi 


N 

■E 

s=l 
1-N-s 




{N-s)l{s-l)l 


IV— «?- — 

X e PI 


'Nt ^Nt 




pIp^ 


N 

E (iv-s)!(s-l)! Vpi 


PI 


-N-s 


„N-S-- 


(80) 


After some algebraic manipulations and with the help of [24, 
Eq. ( 8.432.7)], ( i79b can be computed as ( iSTT i shown on the 
top of the next page. 

To this end, substituting dTSl ) and dSTI ) into dTSb yields the 
desired result. 


Appendix V 
Proof of Theorem[5] 


When pi ^ oo the outage probability of the system can be 
approximated as 

P^out^ ~ Prob (7rT^ < 7th) + Prob (7^2^^ < 7th) ■ (82) 


The c.d.f. of 7^1'^*" can be expressed as 
P^MRc (7ih) = Prob (^llhif^ < , (76) 


where U = 


variables with parameter 


which is an exponential random 
[34]. Hence, we have 


(i-e)p/ 

dZ 


_ ^iT^th AC —1 

dje o-e)pi „ 1 

P,MRC (7th) = 1 - 2^ - 


c^[7th 


poo 

/ 

^0 


1)^ 


■^1 T'th I V 

(i-e)Pi J- 


dx. (77) 


Then, applying the binomial expansion and invoking [24, Eq. 


This approximation comes from the fact that, as pi increases, 
Prob (7^'^'” < 7th and 7^'^*' < 7th) is negligible compared 
with Prob (7^'^^ < 7th) or Prob (7/^'^*' < 7th) (it can be 
justified through Pig |7(b)| i. 

Therefore, the high SNR approximation for the outage 
probability can be given by 


P/out” ~ (7th) + P)^RC (7th) , (83) 

where F°Sirc and F“rc denote the high SNR approximation 

'^11 '^12 

of Prob (7^1'^^ < 7th) and Prob (7)^'^^ < 7th), respectively. 

We start with the characterization of F^rc- Starting from 
( i76b , and with the help of the asymptotic expansion of incom¬ 
plete gamma function, we have 


E)^rc (7th) 


diitu Y^ iii-o)pir 
^jriN-ny.- 


(84) 
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(7th) = 1 - 


A UjJi (1 - ^ - jV rfl dl 

PiPf ^ {N-s)\{s-iy. \pi 


dldljth '' 
vdpi . 

Pi 


Pi 


l-N-s N-1 


N+l-s 


n + s-N~l 



2<-<- ^ UjZl {l-N- j) 
PiP^ h {N-s)\{s-\)\ 


m—Q 


dld}jth \ ' 
■qOpi ) 


n + s-AT-l 


K, 


■,-N-^ 24 


m+s — N—l 


j dldjjth 

TlOpi 


( 81 ) 


Now, we turn our attention to (7th)- According to 

(TT^ and utilizing the the asymptotic expansion of incomplete 

gamma function, conditioned on yi = ^ ||hi||^ and yi = 

^ ||h/||^, F^rc ( 7th) can be expressed as 
^ ’12 


F^rc (7th) 


1 f 47th \ 

r{N + 1) \ye {yi + yi)) 


(85) 


Averaging over yi, we have 


F^rc (7th) 


1 

r{N + i)riN)\ yOpi ) 



1 


yi + x 


N 


dx. 


( 86 ) 


Make a change of variable yi+x = t, and apply the binomial 
expansion, (l86l) can be rewritten by 


e 


( ^2 A 
V ) 


N 


N-1 


f^mrc ( 7 th) r (Af +1) r {N) 


E 


N-1 


i-yi) 


N-i-l 


Pl_ 

d} 

Further averaging over yi, we have 


i-2N+l 


r i-Af + 1, 


d]yi 

PI 


(87) 


poo .4 1 f dldl-ftu 

^MRC th) p p(^)2 


N 


N-1 

E 


N-1 


(- 1 ) 


N-i-l / pi 


J2N-i-2 


i-2N+l 


r ( i - TV + 1, ^ ) dx . (88) 


Xi 


With the help of [24, Eq. (6.455.1)], the integral Fi can be 
solved as 


It = 


/ \ N-i-l / \ N 

(t) (t) r(Jv) 


2W - T - 1 


2Fi(^l,TV;2TV-i;l-^^. (89) 


Then, utilizing [24, Eq. (9.131.1)], we can obtain 


FT]3Rc(7th) 


/ ^1^27th A jv_t 

\ vSpi ) /TV — 1 

r(Tv + i)r(TV) ^ I * , 


(- 1 ) 


AT-i-l 


,Fi (tv, 2W - T - 1; 2TV - T; 1 - gf 
2N-i-l 


(90) 


To this end, substituting ([84]) and ( i90i ) into ( 1831 ) yields the 
desired result. 


Appendix VI 
Proof of Theorem]!] 

Similar to the proof of Theorem]!) we note that the ergodic 
capacity upper bound can be computed as 


/^MRC _ ^ 

^lup — 


^MRC + C^MRC — 
'Tl ’12 


llog,(l + e=('"’"') + eE('"-'“)), ,91) 


where C^urc = [log2 (l + , for fc € {1, 2}. Hence, 

the remaining task is to compute C^mrc, C^URC, E 
and E (in 7^'^'-). 


A. Calculation of C^urc 

■> 7/1 

Utilizing the same method as in [33] and invoking the c.d.f. 
of 7^4'^'" given in (TtHT i. C^uRc can be computed as 


C'.^MRC 

7/1 


1 d\ gg ((l-d)p,f 

2^’^2g7(l-6»)piy ^^dy {m-n)\ 



_ 

g {l-6>)pi 2 ;’^ 


(1 + x) 





I2 


(92) 


With the help of the identity (1 + /3x)~“ = 

pgGig/3x |J““) and the integral formula [36, Eq. 
(2.6.2)], X 2 can be computed as 


12 


((l-0)pi)'"+' 


(i-e)pi 

‘*1 

(l-9)p/ 



■ (93) 


B. Calculation of C„mrc 

Similarly, with the help of the c.d.f. of 7^'^'" given in dSTT) . 
U MRc can be computed as ( l94l i shown on the top of the next 
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(7_mrc = 
112 


pfpf ln2 


N 


rij^i {i-N- j) 

^ (iV-s)!(s- 1)! \pi Pi 

N 


-N-sN-1 ^ 


m—0 


ml 


po) 


AT+l- 


A UUl i^-N-j) Ul d-^j 

pfpfln2^ (iV-s)!(s-l)! \pi pi 


+ 


l-N-s N-1 


E 

ra—Q 


1 

m! 


N + l — S 


l4, (94) 


page, where 


la = 


OO ^N+l-S 


I 


+ X \ pOpi 


pOpi 


and 


r°° rpN+l-s / P'rjT„\ IS:±S—^L-k. 

I4=/ 

Jo 1 + a; V p&pi 


Invoking the c.d.f. of 7/]^^'", the general moment of 


can be computed as 

7V-1 


E((7 


MRC^fe' 


= E 




(1-0) Pi 

j e~T^^x^+^-^n + ^xj dx. 


{m — n)\ 

-(n+l) 


X5 

With the help of [24, Eq. (9.211.4)], we obtain 

7V-1 


E((7,“r)*) = i: 


m—0 


^((1-0^^ 


(l-0)piy {m-n)\ 


k 


d[p7 


m+fc 


r(m + fc)x 


4' m + fc, m + fc — n;' 




(1 - 0)p7 

ed using I 

. To proceed, we find it convenient to use 


fc=o 


It is easy to show that 

dsi (fc) 


dk 


= ln((l - 0) pi) - lnd[ + '0 (1) - 


k^O 


, dTdZx\ 

I 1 dx, (95) 




1,1 
0 , 0,1 I ’ 


2.3 V (1-0)0/ 

where we have utilized the identity T ^1, (i-B)pi 

r{a, 

da 


( 101 ) 


Cl J* 

e , and the derivative property = r (a, 2) In z + 


G5 


(z ^ 

.3 Y 0,0,a J ■ 


( dkdZx\ 

K^+s-n-1 ( 2W J dx. (96) 

Then, following the similar lines as in the Appendix HIH G.^mrc 
can be expressed in closed-form. 

C. Calculation o/E (107^^^^) 

„MRC 


As for 


ds2{k) 

dk 


k=q 

. . dkTi(k) 

to compute — 


dkTi (fc) 


dk 


, it is easy to observe that the key task is 
. Hence, we have 

dTi (k) 


k=0 

= Tiik)\^^Q + k 


k=0 


dk 


( 102 ) 


Noticing that when m > 1, 


MRC 


we have 

Therefore, we get 


fc =0 


= 0, hence 


^ —Q 

dkTi{k) 

dk 


fc=0 

< cx) is a constant, 
fe=o ^ 


ds2 {k) 


dk 


N — 1 m 

= EE 


((1-0)0/)^ 


(97) 


E—f i(n—m)T ( \i ■ 

fe=0 „,=1„=O0/ (m-n)! 

/ J, 

r (to) 4/ I TO, TO — n; ■ 


■ (98) 


The expectation of In 7^'^'- can be derived using E(lnx) = 
dE(a;'“) 


(shown on the top of the next page) as an alternative expression 
of® 

where Ti (fc) = 4- (to + fc, TO + fc - n; (gi) 

xr(m + A:). Now, the expectation of In 7^'^'- can be 
computed as 


MRCt _ (^) 


E (In 7^1'^'^) = 


dk 


ds2 (k) 


k=0 


dk 


( 100 ) 


k=0 


(1 -.)«,'■ 

To this end, substituting (llOlb and (11031) into (llOOl i. the 
expectation of In 7)!'^^ can be obtained. 

D. Calculation of Jl] (In 

The expectation of E (in 774*^*") can be computed as 

E (11177^2'^’“) = lnp0 - IndJ + E (In 2') + E (in ||h2||^^ . 

Since E (in ||h2||p^ = 0 (iV), the remaining task is to figure 
out E(lnZ). With the help of the p.d.f. of Z given in dSOl) . 
we have (1104b shown on the top of the next page. 

Then, invoking [24, Eq. (4.352.1)], the integral in (1104b can 
be solved. 
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